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physical  processes  in  micropropulsion  systems  and  laser-surface  interactions.  The  impulse 
balance  has  been  designed  and  tested  with  a  robust  calibration  system  to  measure  impulsive 
forces  with  resolution  as  low  as  several  nano-Newton-seconds.  The  behavior  of  the  impulse 
balance  was  thoroughly  studied  and  characterized.  A  simple  analytical  model  of  the  balance’s 
motion  was  developed  from  the  general  equation  of  motion  of  an  underdamped,  harmonically 
oscillating  system.  Also,  two  distinct  methods  of  analyzing  the  experimental  data  from  the  nano¬ 
impulse  balance  have  been  investigated.  The  first  method  resolves  the  total  impulse  as  a 
function  of  the  balance’s  maximum  deflection.  The  second  method  enables  the  determination  of 
the  impulse  and/or  force  applied  as  a  function  of  time  from  the  balance’s  time-resolved  motion. 
A  calibration  scheme  employing  electrostatic  actuation  techniques  is  used  to  experimentally 
validate  the  model  and  impulse  measurement  techniques. 
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I.  INTRODUCTION 


Several  impulse  balanees  have  been  designed  to  measure  transient  forees  produeed  by  a  variety 
of  proeesses.^^'^^  Generally,  the  impulse  generated  by  a  deviee  is  investigated  by  measuring  the 
maximum  defleetion  of  the  balanee  eaused  by  a  given  impulse. There  are  two  major 
shorteomings  of  this  approaeh.  First,  there  is  no  information  that  ean  be  inferred  from  this  data 
regarding  the  pulse  width,  r,  or  the  pulse  shape.  Seeond,  it  is  only  strietly  valid  for  t  «  T.  Sinee 
very  different  pulse  widths  and  shapes  ean  lead  to  the  same  maximum  defleetion  on  an  impulse 
balanee,  the  simple  method  of  investigating  a  balanee’s  maximum  defleetion  is  not  adequate  to 
determining  the  eharaeteristies  of  the  impulse  delivered.  A  detailed  analysis  of  the  impulse 
balanee’s  equation  of  motion  ean  lead  to  a  more  eomplete  understanding  of  the  eomplex  physieal 
proeesses  whieh  may  be  driving  deviee  operation. 

There  are  a  number  of  applieations  of  eurrent  interest  where  the  knowledge  of  time  resolved 
forees  would  be  benefieial  to  the  basie  understanding  of  physieal  proeesses.  Two  of  these 
applieations  are  spaeeeraft  propulsion  and  laser  ablation.  For  spaeeeraft  propulsion,  fine  attitude 
eontrol  requires  preeise  and  reprodueible  impulse  delivery  from  the  thruster  system.  Although 
total  delivered  impulse  is  generally  measured  and  deemed  adequate  for  most  mission  seenarios, 
time  resolved  thrust  measurements  would  be  benefieial  in  understanding  the  physies  of  thruster 
operation  and  eould  lead  to  the  development  of  more  preeise  and  efficient  thruster  systems.  For 
example,  pulsed  plasma  thrusters  (PPT)  use  a  spark  generated  across  a  Teflon  propellant  to 
produce  an  energetic  arc-ablated  ionized  plume. The  total  impulse  delivery  from  a  PPT 
involves  a  relatively  low  mass  of  highly  energetic  ions  and  a  large  mass  of  slow  moving 
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molecules  and  heavier  particles  as  shown  schematically  in  Fig.  la.  The  initially  ablated  material 
produces  the  majority  of  the  thrust.  Subsequent  heat  adsorption  by  the  Teflon  produces  a  late¬ 
time  ablative  plume  of  relatively  massive,  slow  moving  molecules  and  particles.  The  late-time 
ablation  does  not  produce  significant  thrust;  however,  a  relatively  large  amount  of  propellant 
mass  is  lost  in  this  process  leading  to  thruster  inefficiencies.  Although  the  total  impulse,  J, 
obtained  from  integrating  F(t)  versus  time  in  Fig.  la  and  Fig.  lb  are  the  same,  the  detailed 
physics  behind  the  production  of  each  impulse  curve  may  be  quite  different.  A  simple  analysis 
based  on  the  balance  maximum  deflection  would  not  yield  any  information  about  the  pulse 
shape.  The  investigation  of  the  time  dependent  forces  produced  by  a  PPT  can  lead  to  more 
efficient  thrusters  by  investigating  system  changes  that  maximize  the  efficient  production  of 
ionized  propellant  and  minimize  the  production  of  late-time  ablation.  A  measure  of  propulsive 
efficiency  is  the  specific  impulse,  fip,  given  by 

(1) 

/=o 

where  F(t)  is  the  time  dependent  force  produced  by  the  thruster,  f  is  the  thruster’s  pulse  duration, 
and  Mp  is  the  total  mass  of  the  propellant  lost  in  the  pulse,  and  go  is  the  Earth’s  gravitational 
constant.  From  Eq.  (1),  it  can  be  seen  that  details  of  the  impulse  delivery  in  the  integral  are 
important  for  thruster  efficiency. 

Within  the  scope  of  laser  ablation,  investigations  of  a  wide  variety  of  physical  processes,  which 
include  sputtering,  vaporization,  ionization,  and  gas  desorption,  are  of  interest. The  conditions 
under  which  each  physical  process  might  occur  and  the  parameters  which  may  control  each 
process  can  be  studied  by  time  resolved  force  measurements.  Eor  example,  an  investigation  of 
the  onset  of  sputtering  as  a  function  of  laser  energy  or  pulse  duration  can  be  performed.  At 
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relatively  low  laser  energy  and  long  pulsewidth,  gas  desorption,  due  to  the  rapid  heating  of  the 
surface,  is  expected  to  be  the  dominant  impulse  producing  process.  As  the  laser  energy  increases 
and/or  the  pulse  duration  decreases,  material  sputtering  or  ionization  may  become  dominant.'-^^ 
Through  a  careful  investigation  of  the  time  resolved  forces  produced  in  the  photon-surface 
interactions,  the  energy  corresponding  to  the  onset  of  sputtering  for  a  particular  material  may  be 
identified. 

The  scope  of  this  work  was  three  fold.  First,  a  uniquely  successful  nano-Newton-second  impulse 
balance  system  (NIBS)  capable  of  resolving  total  impulse  measurements  as  low  as  7  nano- 
Newton-seconds  (nNs)  was  designed  and  constructed.  Second,  an  analytical  model  was 
developed  from  the  equations  of  motion  for  an  underdamped,  harmonically  oscillating  system. 
The  model  was  validated  using  experimental  data  obtained  from  the  NIBS.  Finally,  a  data 
reduction  method  to  accurately  time-resolve  the  applied  force  from  the  NIBS  deflection  data  was 
developed.  The  developed  model  and  data  reduction  methods  in  this  study  are  extended  to 
impulse  balance  operation  for  arbitrary  pulsewidth  (i.e.  with  no  restriction  on  the  pulsewidth 
duration  relative  to  the  natural  motion  of  the  stand). 

II.  NANO-IMPULSE  BALANCE  SYSTEM 


The  current  version  of  the  NIBS  is  based  upon  a  torsional  thrust  stand  developed  jointly  by  the 
Air  Eorce  Research  Eaboratory  and  the  University  of  Southern  California. The  NIBS  design 
has  been  modified  to  allow  for  very  low  impulse  measurements,  which  could  not  be  attained  by 
previous  versions  of  the  stand.  An  electrostatic  comb  force  calibration  technique  described  by 
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Selden  and  Ketsdever^^^  has  been  employed  to  aeeurately  calibrate  the  system  for  transient 
impulses. 

As  shown  in  Fig.  2,  the  NIBS  is  a  torsional  pendulum  with  viscous  damping.  Two  flexure  pivots 
provide  the  restoring  force  for  the  system  with  a  combined  rotational  spring  constant  of 
approximately  0.200Nm/rad  (measured).  The  entire  structure  is  made  of  aluminum  to  keep  the 
balance’s  moment  of  inertia  as  low  as  possible  in  an  attempt  to  maximize  the  deflection  for  a 
given  applied  force.  The  arms  were  fabricated  from  square  tubing  to  provide  rigidity  over  long 
lengths,  and  channels  were  cut  from  the  side  of  the  tubing  to  reduce  mass.  Propulsion  systems  or 
test  materials  are  mounted  at  the  end  of  the  46  cm  long  arms.  The  Electrostatics  Force 
Calibration  System  (EFCS)  is  mounted  on  the  end  of  one  arm  at  the  same  location  from  the 
center  of  rotation  as  the  impulse  delivery  is  expected.  The  NIBS  motion  is  sensed  by  a  linear 
variable  differential  transducer  (EVDT).  The  EVDT  senses  a  linear  motion  at  the  end  of  the 
balance’s  arm,  which  can  be  transformed  into  a  rotational  deflection.  The  data  trace  obtained  by 
monitoring  the  EVDT  output  over  time  provides  the  basis  for  deriving  the  impulse  imparted  to 
the  NIBS.  A  damping  cup  is  mounted  to  the  bottom  of  the  NIBS  and  is  placed  in  an  oil  bath  to 
provide  viscous  damping,  which  is  essential  for  damping  out  the  motion  of  the  stand  in  vacuum. 

III.  IMPULSE  BALANCE  MODEL 


The  motion  of  an  underdamped,  harmonically  oscillating  system  can  be  expressed  in  rotational 
or  linear  terms  by  the  following 

Wit)  +  Ceit)  +  KBit)  =  Fit)r  =  Mit)  (2) 
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Ix{t)  +  Cx{t)  +  Kx{t)  =  F{t)r  ^  =M{t)r 


The  angular  solutions  can  easily  be  transformed  into  a  linear  deflection  to  correspond  to  the 
measurement  methods  used  in  the  experimental  setup.  For  the  general  case  of  forced  motion  of 
an  underdamped  system,  consider  M(t)  =  Mo  to  be  constant  over  a  defined  time  period  (pulse 
width)  T,  as  shown  in  Fig.  lb.  This  can  be  representative  of  either  a  constant  impulse  or  steady- 
state  situation,  depending  on  the  length  of  t.  The  solution  to  Eq.  (2)  for  initial  conditions  of 
angular  displacement  Oo  and  angular  velocity  <9^ 
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where  a  =  -C/21  =  -5/T  and  P  =  [K/I  -  a  ]  .  Note  that  by  removing  the  damping  term,  the 

1  /2 

expected  natural  frequency  cOn  =  [K/I]  for  a  simple  torsional  spring  is  obtained. 


A  closed-form  analytical  solution  of  the  equation  of  motion  in  Eq.  (2)  is  not  possible  for  any 
arbitrary  applied  moment,  M(t).  However,  the  problem  may  be  broken  up  into  segments.  Any 
arbitrary  or  irregular  force  can  be  reasonably  approximated  using  a  series  of  constant  force 
segments  of  small  time  widths  dr.  The  initial  conditions  of  any  given  segment  are  simply  the 
end  state  of  the  displacement  and  velocity  terms  from  the  previous  segment. 
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Equations  (3-5)  have  been  used  to  develop  a  simple  software  tool  that  has  been  validated  using 
experimental  data.  The  software  has  been  used  to  examine  the  expected  behavior  of  the  NIBS 
for  impulses  of  varying  pulse  widths  and  for  steady  state  operations.  Figure  3  shows  a 
comparison  of  NIBS  experimental  data  for  a  242  nNs  impulse  with  the  analytical  model.  Figure 
3  shows  good  agreement  between  the  model  and  experimental  data  for  the  magnitude  and 
frequency  of  the  balance’s  motion. 

IV.  SIGNAF  PROCESSING  AND  DATA  REDUCTION 


A.  Wavelet  And  Fourier  Transform  Denoising 

The  ability  to  obtain  reliable  and  repeatable  force  measurements  from  NIBS  depends  on  the 
time-varying  deflection  traces  provided  by  the  EVDT.  The  EVDT  sensor  provides  an  analog 
output,  which  commonly  exhibits  significant  high  frequency  noise  components.  The  noise  can 
mostly  be  attributed  to  electrical  sources  (primarily  in  the  15Hz,  50Hz  and  60Hz  ranges).  Two 
signal  analysis  techniques  have  been  employed  in  this  study  to  perform  denoising  of  the  EVDT 
signal:  Fourier  and  wavelet  transforms. 

First,  the  Fourier  transform  F(co)  of  a  given  signal  f(t)  deconstructs  the  signal  into  constituent 
sinusoidal  components.  Fourier  signal  analysis  techniques  are  fairly  common  and  numerous  tools 
exist  to  facilitate  the  transform  of  a  given  signal  into  a  frequency  spectrum.  The  denoising 
technique  employed  was  a  simple  frequency  thresholding,  whereby  all  components  higher  than  a 
specified  frequency  are  filtered  out.  This  technique  essentially  functions  as  an  ideal  low-pass 
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filter  (LPF)J^'^  The  signal  can  then  be  transformed  back  from  the  frequency  domain  into  the 
time  domain  for  further  data  reduction  and  analysis. 

While  the  principal  frequency  of  interest  is  the  effective  frequency  of  the  NIBS,  care  must  be 
taken  in  choosing  a  threshold  frequency.  Selecting  too  low  of  a  threshold  frequency  may  result 
in  the  loss  of  actual  signal,  which  would  typically  result  in  the  undervaluing  of  the  derived 
impulse.  For  sets  of  data  that  were  analyzed  using  the  Fourier  techniques,  the  threshold 
frequency  used  was  10  Hz. 

Alternately,  the  wavelet  transform  has  also  proved  to  be  a  valuable  tool  for  denoising  and 
smoothing  of  the  raw  experiment  data.  The  wavelet  transform  is  analogous  to  Fourier  transforms 
in  that  it  decomposes  a  signal  into  coefficients  of  constituent  waveforms.  However,  Fourier 
analysis  depends  entirely  on  decomposition  by  sinusoids  across  the  frequency  spectrum. 
Wavelet  analysis  allows  for  the  waveform  to  be  of  any  arbitrary  shape  and  decomposes  the 
signal  across  the  time  domain.  The  coefficients  of  the  wavelet  transform  are  a  function  of  scale 
and  position  (shift  in  time).^^^^ 

The  decomposition  by  scale  and  position  allows  for  the  isolation  of  successive  levels  of  small 
scale  structure  or  detail,  leaving  the  remaining  portion  as  an  approximation  of  the  original  signal. 
Applying  thresholding  settings  to  these  detail  coefficients  can  result  in  removing  general  noise 
components  while  leaving  actual  fine  structure  intact.  Such  a  deconstruction  into  scale  and  time 
domains  can  provide  valuable  information  about  the  signal  through  inspection  of  the  coefficients, 
as  anomalous  points  and  localized  events  are  often  clearly  evident.  Also,  since  the  transform 
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remains  in  the  time  domain,  it  is  also  possible  to  perform  wavelet  analysis  on  finite  or  shifting 
segments  of  a  data  set.  This  can  be  an  extremely  useful  in  handling  large  data  sets.  The  major 
drawback  of  the  wavelet  process  is  that  it  can  be  extremely  computationally  intensive. 

Both  transform  methods  have  been  demonstrated  to  be  quite  reliable  for  deriving  the  total 
impulse  from  the  deflection  data.  In  either  case,  attention  must  be  paid  to  not  filter  out  more  than 
just  the  noise  from  the  signals.  Figure  4  shows  a  comparison  between  the  Fourier  and  wavelet 
denoising  techniques  for  a  typical  deflection  trace.  Agreement  between  the  Fourier  and  wavelet 
processed  signals  is  excellent.  Furthermore,  the  comparison  of  the  derived  impulse  from  such 
traces  processed  with  Fourier  and  wavelet  techniques  resulted  in  differences  less  than  1%. 

B.  Time-Resolved  Impulse  Measurement  Process 

Once  the  raw  data  has  been  denoised  through  either  the  Fourier  and/or  wavelet  transform 
methods,  it  can  be  used  to  derive  the  total  impulse  that  was  imparted  to  the  NIBS.  In  the 
simplest  case,  for  r  «  T,  the  total  impulse  may  be  derived  from  the  maximum  deflection. 
Calibration  curves  for  this  method  can  readily  be  derived,  as  will  be  shown  in  Section  VI. 
However,  the  real  interesting  details  of  the  impulse  or  force  applied  are  the  time  resolved 
characteristics.  The  Time -Resolved  Impulse  Measurement  (TRIM)  process  is  derived  from  Eq. 
(2).  The  equation  of  motion  for  the  NIBS  indicates  that  the  force  applied  (as  a  function  of  time) 
is  simply  related  to  the  addition  of  the  position,  velocity  and  acceleration  components,  which  are 
scaled  by  the  spring  constant  K,  damping  coefficient  C,  and  the  moment  of  inertia  I,  respectively 
as  shown  in  Fig.  5. 


9 


Since  the  position  (deflection)  as  a  function  of  time  is  measured  in  the  experiments,  the 
derivative  of  the  data  must  simply  be  taken  onee  to  obtain  the  velocity  and  twice  for  the 
acceleration.  The  major  problem  here  is  resolving  the  derivatives  for  a  signal  with  noise.  Using 
the  signal  processing  methods  previously  described,  the  raw  position  data  can  be  cleaned  up 
significantly.  In  addition,  a  Savitzky-Golay  algorithm  is  used  to  take  the  derivatives.  This 
algorithm  employs  a  moving  window  over  whieh  the  points  are  fitted  for  a  line.  The  slope  of 
the  line  represents  the  derivative  of  the  point  at  the  center  of  the  window.  This  method  results  in 
artifacts  at  the  beginning  and  end  of  a  data  trace.  This  method  also  tends  to  introduce  slight 
aliasing  effects,  depending  on  the  size  of  the  window.  However,  it  reliably  maintains  the 
content  under  the  curve,  which  is  essential  for  ealculating  the  total  impulse. 

To  complete  the  proeess,  we  need  to  determine  I,  C  and  K.  The  moment  of  inertia  of  the  NIBS 
can  be  determined  through  a  combination  of  physical  measurement  and  CAD  analysis.  The 
spring  eonstant  and  damping  coefficient  ean  then  be  determined  from  the  position  data  by 
extraeting  the  period  of  the  oscillations  T  and  the  damping  ratio  5,  which  is  obtained  from  the 
natural  logarithm  of  the  ratio  of  successive  peaks  in  the  displacement  data,  5  =  ln[xmaxi/xmax2]. 
For  completeness,  it  is  also  necessary  to  remember  to  include  the  term  for  the  moment  arm  in  Eq. 
(2).  Finally,  from  the  derived  force  F(t),  the  total  impulse  J(t)  is  given  by  integrating  F(t). 

V.  EXPERIMENTAF  SETUP  AND  PROCEDURE 


The  NIBS  was  installed  in  a  41  cm  diameter  x  122  cm  long,  stainless  steel  vaeuum  ehamber 
fitted  with  a  450  E/see  turbomoleeular  pump  capable  of  maintaining  pressures  of  approximately 
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2  X  10'^  Torr.  For  this  study,  impulse  delivery  to  the  NIBS  is  aeeomplished  by  supplying  a 
potential  difference  to  the  EFCS.  The  EFCS  consisted  of  either  an  aluminum  parallel  plate 
assembly'-^^  or  an  aluminum  comb  assembly'^^^  with  one  side  of  the  EECS  assembly  attached 
directly  to  the  NIBS.  As  shown  in  the  experimental  setup  in  Eig.  6,  the  power  supply  for  the 
EECS  was  attached  to  a  pulse  generator  capable  of  delivering  ±3500  V  with  a  20  ns  rise  time  and 
a  variable  DC  pulse  width  (minimum  of  60  ns).  The  output  of  the  pulse  generator  was  sent 
directly  to  the  EECS  assembly  in  the  chamber  through  a  high  voltage  vacuum  feedthrough  and 
was  monitored  by  a  16-bit  data  acquisition  (DAQ)  system  with  a  maximum  sample  rate  of 
333kHz.  The  applied  (and  monitored)  voltage  to  the  EECS  was  used  to  calculate  the  actual  time- 
dependent  force  applied  to  the  NIBS.'-^^ 

The  actual  forces  applied  were  later  compared  to  the  force  derived  from  the  NIBS  deflection  data 
through  the  TRIM  process  described  in  Section  IV-B.  The  motion  of  the  NIBS  was  measured 
using  an  EVDT,  which  was  connected  to  either  a  16  or  24-bit  data  acquisition  system.  The  24- 
bit  system  used  was  limited  to  a  sampling  rate  of  60  Hz,  and  the  16-bit  system  was  capable  of 
sampling  at  333  kHz.  In  some  of  the  tests,  a  simple  low-pass  filter  (EPE)  was  used  on  the  EVDT 
signal  to  minimize  electrical  noise.  The  analog  filtering  serves  to  complement  the  signal 
processing  techniques  discussed  in  Section  IV-A.  The  raw  data  was  fdtered  to  reduce  the  noise 
in  the  experimental  data,  and  the  denoised  data  was  differentiated  to  provide  NIBS  velocity  and 
acceleration  information,  which  was  used  to  determine  the  characteristics  of  the  delivered 
impulse  as  described  in  the  previous  section.  Values  of  T  and  5  were  obtained  using  NIBS 
displacement  data.  Erom  these  parameters  and  the  measured  moment  of  inertia,  the  values  of  the 
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damping  (C)  and  spring  (K)  coefficients  were  calculated.  Table  1  shows  the  characteristic  values 
for  the  NIBS  in  this  study. 

VI.  RESULTS 

A.  Calibration 

Calibration  of  the  NIBS  was  a  critical  element  of  this  study.  The  application  of  a  variety  of 
known  impulses  of  different  magnitudes  and  pulsewidths  were  tested  to  demonstrate  the  stability 
and  repeatability  of  the  NIBS.  From  these  tests,  a  process  has  been  developed  which  allows  the 
balance  to  be  well  characterized  under  any  configuration. 

In  order  to  use  the  LVDT  voltage  output  to  obtain  a  derived  impulse,  a  calibration  factor  is 
required  that  converts  the  analog  voltage  signals  into  appropriate  displacement  measurements. 
The  LVDT  was  independently  calibrated  from  repeated  measurements  in  a  staged  micrometer 
configuration  that  spanned  the  full  range  of  the  LVDT.  The  conversion  factor  was  determined 
to  be  1453.366  V/m  for  the  tested  configuration,  which  included  the  use  of  a  LPF  for  analog 
signal  conditioning.  For  test  configurations  without  the  use  of  a  LPF,  an  additional  scaling  factor 
of  1.0902  was  necessary  to  compensate  for  the  signal  attenuation  caused  by  the  LPF. 

The  applied  impulse  is  the  known  impulse  delivered  to  the  stand  using  the  LEGS.  The  voltage 
applied  to  the  EFCS  is  typically  greater  than  lOOV  and  is  monitored  and  known  to  within  O.IV. 
The  actual  force  produced  by  the  EFCS  was  independently  calibrated  by  the  microbalance 
method  described  by  Selden  and  Ketsdever.^^^  The  error  for  a  given  applied  force  was  found  to 
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be  less  than  2%.  For  all  results  presented  in  this  paper,  this  error  is  understood  to  be  present  in 
all  applied  forces  and  impulses,  whether  or  not  it  is  visible  in  the  figures  or  explicitly  stated  in 
the  values. 

Figure  7  shows  three  sets  of  dissimilar  calibration  data  that  was  analyzed  using  appropriate 
calibration  factors.  One  set  of  data  was  originally  sampled  at  9000Hz,  the  next  was  lOOOHz 
data,  and  the  third  was  lOOOHz  with  the  use  of  an  analog  LPF.  Each  of  the  sets  show  reasonable 
agreement  with  expected  results  as  the  derived  impulses  are  consistently  within  4.5%  of  the 
applied  impulse.  Each  calibration  data  set  also  shows  trends  of  being  consistently  slightly  high 
or  low,  which  would  be  consistent  with  filtering  variations  for  each  of  the  sets  of  data.  This 
suggests  that  for  any  configuration  of  the  NIBS,  a  consistent  set  of  calibration  data  could  be  used 
with  a  consistent  signal  processing  parameter  to  obtain  appropriate  calibration  factors  for  an 
unknown  impulse,  which  could  minimize  error  in  the  derived  impulse.  Nevertheless,  without 
any  additional  correction  factors,  the  results  for  derived  impulse  throughout  the  paper  repeatedly 
fall  within  4.5%  of  the  applied  impulse. 

B.  Maximum  Deflection  Analysis 

Eigure  8  shows  the  NIBS  maximum  linear  deflection  as  a  function  of  applied  impulse.  The 
analytical  model  is  in  good  agreement  with  the  experimental  data  as  shown  in  Eig.  8.  Eor  the 
cases  presented  in  Eig.  8,  the  impulse  pulsewidths  (t)  were  shorter  than  approximately  one-tenth 
of  the  period  (T)  of  the  NIBS.  Typically,  impulse  balances  are  operated  using  the  assumption 
that  the  maximum  deflection  is  linearly  dependent  on  the  delivered  impulse. The  assumption 
is  only  strictly  valid  for  r  «  T.  Eigure  9  shows  the  maximum  NIBS  deflection  as  a  function  of  r  / 
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T  for  various  total  impulses.  The  analytieal  model  eaptures  not  only  the  trend  of  the 
experimental  data,  but  is  also  in  excellent  agreement  with  the  magnitude  of  the  data.  As  seen  in 
Fig.  9  for  T  >  0.1  T,  the  direct  correlation  between  the  maximum  balance  displacement  and  the 
delivered  impulse  is  no  longer  valid.  For  r  «  T,  the  maximum  deflection  has  a  linear  relationship 
with  the  total  impulse  as  shown  in  Fig.  8.  For  larger  t,  there  are  two  possible  scenarios.  First, 
two  different  total  impulses  (say  one  with  r  «  T  and  one  with  r  ~  T)  can  lead  to  the  same 
maximum  deflection;  second,  two  identical  total  impulses  with  different  pulsewidth  can  result  in 
different  maximum  deflection.  Examples  of  these  two  scenarios  will  be  shown  in  Section  VI-C. 
Therefore,  another  means  by  which  to  analyze  the  data  is  necessary  if  the  pulsewidth  is  not 
known  a  priori. 

C.  Derived  Impulse  Analysis 

Figure  10a  shows  the  impulse  balance  deflection  as  a  function  of  time  for  two  impulses  of  the 
same  total  magnitude  (integrated  force  with  time).  The  maximum  deflection  from  the  NIBS  is 
not  the  same  even  though  the  total  impulse  is  identical  for  both  cases.  The  actual  time  varying 
force  applied  is  shown  by  the  dashed  line  in  Fig.  10b.  As  shown,  the  pulsewidth  of  the  two 
traces  was  different  by  an  order  of  magnitude.  The  TRIM  process  described  in  Section  IV-B  was 
used  to  obtain  the  derived  time  varying  force  also  shown  in  Fig.  10b.  Figure  10c  shows  the  total 
impulse  delivered  to  the  balance  obtained  by  the  integral  under  the  curves  in  Fig.  10b.  Although 
the  maximum  deflection  for  the  two  pulses  was  not  the  same.  Fig.  10c  shows  that  the  total 
impulse  delivered  were  in  fact  identical. 
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Figure  11a  shows  the  NIBS  linear  displacement  versus  time  for  impulses  of  42.34  pNs  and  21.65 
pNs.  The  pulsewidths  for  these  impulses  was  1.55  sec  and  0.1  sec  respectively.  As  seen  in  Fig. 
10a,  the  maximum  NIBS  deflection  is  the  same  for  both  impulses  even  though  the  total  impulses 
are  different  by  a  factor  of  two.  Figure  10b  shows  the  impulse  shape  for  both  cases.  The  dashed 
lines  are  from  the  signal  delivered  by  the  pulse  generator,  and  the  solid  lines  are  derived  from  the 
filtered  NIBS  deflection  data  using  the  TRIM  data  reduction  technique  described  in  Section  IV- 
B.  As  Fig.  11b  indicates,  there  are  system  time  response  issues  that  would  not  allow  better 
resolution  of  the  impulse  shape  with  faster  data  acquisition  sampling.  The  time  response  of  the 
LVDT  signal-conditioning  unit  was  seen  to  be  less  than  100  Hz,  which  was  generally  not  fast 
enough  to  adequately  resolve  the  impulse  shape  for  very  short  pulsewidths  due  to  aliasing. 
However,  the  total  impulse  taken  as  the  integral  under  both  curves  (pulse  generator,  LVDT 
derived)  is  the  same  to  within  4.5%  as  shown  in  Fig.  lie.  Although  the  maximum  displacement 
of  the  NIBS  was  the  same  for  the  two  impulses  (as  shown  in  Fig.  11a),  the  total  impulse  was 
indeed  a  factor  of  two  different  as  shown  in  Fig.  lie. 

Figure  12a  shows  the  experimentally  derived  force  from  the  NIBS  as  a  function  of  time  for  an 
impulse  comprised  of  two  square  waves  from  the  pulse  generator.  The  pulses  each  represent  a 
force  of  170.7  pN  over  a  1.0  sec  pulsewidth  and  are  separated  by  2.0  sec.  As  seen  in  Fig.  12a, 
there  is  good  agreement  between  the  measured  impulses  delivered  by  the  pulse  generator  to  the 
EFCS  and  the  derived  impulses  from  the  NIBS  deflection  data.  Figure  12b  shows  the  delivered 
impulse  as  a  function  of  time.  As  Fig.  12  shows,  the  data  analysis  method  developed  in  this 
study  is  adequate  to  resolve  discrete  forces  applied  to  the  impulse  balance. 
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Although  attention  has  been  paid  to  only  square-wave  impulses  to  this  point,  arbitrary  impulse 
shapes  ean  also  be  investigated.  Figure  13a  shows  a  time  varying  force  applied  to  the  thrust 
stand  which  might  resemble  that  produced  by  a  PPT  or  laser  ablation  mechanism.  The  derived 
force  obtained  from  the  NIBS  position  versus  time  data  reproduced  the  actual  applied  force 
reasonably  well.  Many  of  the  differences  in  the  time  varying  force  have  been  linked  to 
limitations  in  the  time  resolution  of  the  LVDT  conditioning  unit.  Although  a  factor  of  perhaps 
10  improvement  in  the  LVDT  can  be  reasonably  expected,  other  methods  of  time  resolving  the 
NIBS  deflection  are  being  investigated.  As  shown  in  Fig.  13b,  the  applied  and  derived  impulses 
are  nearly  identical  and  compare  to  within  1.9%. 

D.  Model  Comparison 

Figure  14  shows  the  comparison  of  the  developed  model  with  experimental  data  for  a  relatively 
complicated  case  of  four  square-wave  impulses  applied  over  a  25  second  time  period.  In  Fig. 
I4a,  the  model  is  able  to  capture  the  motion  of  the  NIBS  for  a  series  of  impulses  delivered  to  the 
balance.  The  model  used  the  measured  and  calculated  characteristics  of  the  NIBS  given  in  Table 
1  as  input  parameters.  Figures  14b  and  c  show  the  calculated  values  of  the  NIBS  velocity  and 
acceleration  as  a  function  of  time.  Again,  the  model  compares  well  with  the  experimentally 
derived  data  obtained  from  differentiating  the  NIBS  time  varying  position  data.  Figure  I4d 
shows  the  actual  time  varying  forces  applied  to  the  NIBS  compared  to  the  experimentally 
derived  data.  The  first  pulse  was  98.39  pN  with  a  0.2  sec  pulsewidth  with  the  second  pulse  of 
the  same  magnitude  separated  by  0.3  sec.  The  second  and  third  pulses  were  separated  by  8.606 
sec.  The  third  and  fourth  pulses  were  also  98.39  pN  with  a  0.4  sec  pulsewidth  separated  by 
7.055  sec.  The  integrals  under  the  time  varying  force  curves  are  shown  in  Fig.  14e.  Again,  the 
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experimental  data  was  able  to  time  resolve  the  impulse  delivery  to  the  NIBS,  and  the  total 
impulse  derived  from  the  NIBS  data  is  within  2.2%  of  the  actual  applied  impulse. 


VII.  DISCUSSION 


The  extension  of  operating  conditions  of  impulse  balances  to  arbitrary  impulse  pulsewidth, 
shape,  and  magnitude  (down  to  7  nNs)  has  been  demonstrated  with  the  NIBS  using  a 
combination  of  an  analytical  code  and  a  data  analysis  method.  Previous  studies'-^’^^  have  been 
limited  to  data  analysis  based  on  the  assumption  that  r  «  T.  For  basic  research  involving  new 
propulsion  systems  and  photon-surface  interactions,  the  characteristic  time-dependence  of  the 
impulse  of  a  particular  process  may  not  be  known  a  priori,  making  a  generalized  scheme 
necessary  to  time-resolve  the  impulse  delivery  to  the  balance.  Also,  a  torsional  impulse  balance 
can  deflect  the  same  amount  for  different  total  impulses  (e.g.  see  Fig.  9),  particularly  where  the 
assumption  of  r  «  T  is  no  longer  valid.  Under  similar  pulsewidth  conditions,  an  impulse  balance 
can  also  deflect  differently  for  equivalent  total  impulses.  Therefore,  a  generalized  r-independent 
data  analysis  tool  is  required. 

The  NIBS  has  already  demonstrated  that  it  can  be  used  to  time -resolve  impulses.  The  real  issues 
now  become  that  of  performance.  The  common  measure  of  an  impulse  balance’s  performance  is 
typically  that  of  impulse  magnitude  resolution.  The  NIBS  has  repeatably  measured  impulses  as 
low  as  7nNs  and  as  high  hundreds  of  micro-Newton-seconds,  and  has  consistently  derived  results 
that  fall  within  4.5%  (worst  case)  of  the  applied  impulse.  However,  when  time  resolving  the 
forces,  it  becomes  apparent  that  there  are  more  performance  related  issues  than  just  the 
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magnitude  of  the  impulse.  Qualitatively,  the  aecuracy  of  the  shape  of  the  time-varying  force  is 
of  interest.  In  other  words,  how  closely  the  derived  force  matches  the  applied  force  over  time  is 
important  in  determining  the  overall  effectiveness  of  both  the  TRIM  technique  and  the  NIBS 
itself.  Noise  and  time  response  of  the  whole  system  then  become  significant  issues. 

Noise  in  the  measured  signal  from  the  NIBS  is  one  of  the  largest  sources  of  error  in  the  TRIM 
process.  Excessive  noise  affects  the  derivatives  of  the  deflection  trace  and  results  in  noisy  time- 
resolved  forces  or  impulses.  There  are  several  sources  of  noise  including  external  physical 
vibrations,  temperature  variations,  and  electronic  noise  directly  from  the  LVDT  and  its  signal 
conditioner.  The  simplest  way  to  eliminate  a  significant  portion  of  the  electronic  noise  is  the 
use  of  a  well-tuned  EPF.  However,  EPFs  introduce  significant  time  response  problems,  as  they 
attenuate  high  frequency  components  that  may  actually  be  part  of  the  signal.  Only  when  time¬ 
resolving  the  forces  is  not  necessary,  the  EPF  provides  the  best  method  to  obtaining  clean  signals 
for  basic  maximum  deflection  measurements.  In  all  other  cases,  the  EPF  should  not  be  utilized, 
and  the  signal  processing  discussed  in  Section  IV  becomes  the  critical  noise  reduction  method. 

In  terms  of  system  time  response,  each  component  of  the  NIBS  system  must  be  examined.  First, 
the  DAQs  components  play  an  important  role.  Often  there  is  a  trade-off  between  resolution  and 
speed.  The  NIBS  system  is  equipped  with  both  a  16-bit  resolution  system  that  samples  up  to 
333KSamples/sec  and  a  24-bit  resolution  system  that  samples  up  to  60  Samples/sec.  Ideally,  a 
high-speed,  24-bit  resolution  system  would  be  optimal,  but  has  yet  to  be  identified. 
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Next,  the  impulse  balance  itself  must  be  considered  for  time  response  issues.  The  NIBS  has  been 
extensively  tested  with  impulses  of  pulsewidth  r  as  low  as  10'^  sec  up  to  minutes  in  length.  Due 
to  time  response  limitations  from  the  LVDT,  the  TRIM  process  only  currently  provides  accurate 
results  for  pulsewidths  as  low  as  0.05  seconds.  Below  this  point,  aliasing  effects  make  it  difficult 
to  discern  whether  the  time -resolution  problems  are  due  to  the  impulse  balance  or  just  signal 
processing  artifacts,  though  the  latter  is  suspected  to  be  the  most  likely.  However,  from 
extensive  maximum  deflection  results,  the  NIBS  behaves  as  expected  across  this  entire  range  of 
pulsewidths.  This  suggests  that  the  stand  is  behaving  well  and  does  not  have  time  response 
problems,  at  least  on  the  order  of  10'^  seconds. 

The  last  system  component  to  examine  for  time  response  issues  is  the  LVDT  and  its  signal 
conditioner.  The  time  response  of  an  LVDT  is  typically  considered  to  be  at  best  one  tenth  of  the 
excitation  frequency  of  the  LVDT.  For  the  NIBS,  the  LVDT  is  driven  at  an  excitation  frequency 
of  2.5KHz.  However,  the  signal  conditioning  provided  within  the  controller  unit  limits  the 
response  characteristics  further. This  accounts  for  the  current  time-resolution  limitations  of 
the  TRIM  process.  Improvement  of  the  technique  described  here  would  be  easily  achieved 
through  improvements  of  the  overall  system  time  response.  While  providing  excellent 
displacement  resolution,  the  LVDT  signal  conditioner  was  found  to  have  the  slowest  time 
response  of  the  entire  system,  making  its  replacement  desirable.  Optical  interferometry  may 
offer  improved  time  response  over  the  LVDT  system  and  has  been  previously  investigated, 
although  not  for  these  particular  impulse  derivation  methods.  Alternatively,  capacitive  and  other 
electromagnetic  sensors  could  also  be  considered. 
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For  actual  operation  of  the  NIBS,  the  error  in  the  TRIM  proeess  ean  be  redueed  signifieantly  by 
using  a  consistent  ealibration  scheme.  Sinee  any  modifieation  to  the  impulse  balance  ean  change 
its  fundamental  characteristies  (i.e.  I,  C,  and  K),  it  is  neeessary  to  ealibrate  in  the  aetual  test 
eonfiguration.  Using  a  robust  set  of  ealibration  tests  that  maintains  a  set  sampling  frequeney  and 
eovers  a  wide  range  of  impulses,  forces  and  pulsewidth  variations,  a  eorreetion  faetor  eould  be 
determined  that  would  minimize  the  error.  Ideally,  ealibration  sets  will  be  taken  before  and  after 
any  experimental  measurements  to  ensure  that  the  impulse  balanee  remains  eonsistent. 

Ultimately,  the  performanee  of  the  impulse  balanee  comes  down  to  repeatability.  For  the 
maximum  deflection  analysis,  the  experimental  error  was  taken  as  the  standard  deviation  of 
between  five  and  ten  data  points  and  was  typieally  less  than  ±2%  (i.e.  smaller  than  the  symbol 
size  in  Figs.  8  and  9  unless  otherwise  indieated).  For  the  time-resolved  impulse  analysis,  all 
derivations  of  total  impulse  from  the  deflection  measurements  have  resulted  in  errors  less  than 
±4.5%,  whieh  encapsulates  all  sourees  of  error  ineluding  from  the  applied  foree.  The 
repeatability  of  the  NIBS  data  suggests  that  it  is  a  viable  diagnostic  tool  for  the  investigation  of 
arbitrarily  produeed  impulses  from  a  variety  of  deviees. 

NOMENCLATURE 


a  damping  term  in  equations  of  motion  (see"') 

P  frequeney  term  in  equations  of  motion  (see"') 

5  damping  ratio,  In[xmaxi/Xmax2] 

3  impulse,  j  F {t)dt  (N-sqc) 
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xor  V 
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pulsewidth  (sec) 

angular  displacement  or  deflection  (rad) 
angular  velocity  (rad/sec) 
angular  acceleration  (rad/sec  ) 
effective  (angular)  frequency  (sec'') 
natural  (angular)  frequency  (sec'') 
viscous  damping  coefficient  (kg-m  /sec) 
force  (N) 

moment  of  inertia  (kg-m  ) 
spring  constant  (N/m) 
moment  or  torque,  F  •  r  (N-m) 
moment  arm  (m) 
period  (sec) 
time  (sec) 

linear  displacement  or  deflection  (m) 
linear  velocity  (m/sec) 
linear  acceleration  (m/sec  ) 
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Figure  1 


Figure  2 


Figure  3 


Figure  4 


Figure  5 


Figure  6 


General  NIBS  Charaeteristies 

Sehematie  of  two  very  different  foree  profiles  that  result  in  equivalent  total 
impulses  Ji=J2-  (a)  Farge  initial  foree  with  smaller  residuals,  (b)  Constant  foree 
of  pulsewidth  r. 

Sehematie  of  the  nano-Newton  impulse  balanee  system. 

Time-varying  defleetion  traee  for  a  242  nNs  impulse  eomparing  experimental  data 
and  the  analytieal  model. 

Denoising  of  FVDT  experimental  data  using  wavelet  transform  and  Fourier 
transform  teehniques. 

Sehematie  representation  of  the  TRIM  proeess  for  deriving  the  time-varying 
impulse  or  foree  from  the  balanee’s  defleetion  and  its  derivatives. 

Sehematie  of  NIBS  experimental  setup  and  TRIM  proeess  for  obtaining  the 
derived  and  applied  forees  imparted  to  the  impulse  balanee  as  a  funetion  of  time. 
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Figure  7 


Figure  8 


Figure  9 


Figure  10 


Figure  1 1 


TRIM  derived  impulse  versus  applied  impulse  for  three  dissimilar  sets  of 
ealibration  data. 

Comparison  of  maximum  deflection  versus  applied  impulse  of  the  analytical 
model  and  experimental  deflection  data.  Error  bars  represent  the  standard 
deviation  of  at  least  10  repeated  tests,  and  are  typically  smaller  than  the  symbol 
used. 

NIBS  maximum  deflection  as  a  function  of  normalized  impulse  pulsewidth 
compared  to  the  analytical  model  results.  For  r  «  T,  the  maximum  deflection  is 
constant  .  The  data  points  represent  experimental  measurements  with  error  bars 
showing  the  standard  deviation  of  at  least  10  repeated  tests. 

Two  test  cases  of  equivalent  total  impulse,  with  different  applied  force  and 
pulsewidth.  (a)  shows  that  equivalent  total  impulses  can  result  in  different 
maximum  displacements,  (b)  shows  the  shape  of  the  applied  force  for  each 
impulse,  (c)  shows  that  the  two  total  impulses  are  in  fact  equivalent. 

Two  test  cases  of  different  total  impulses  that  result  in  the  same  maximum 
deflection  of  the  balance  as  shown  in  (a),  (b)  shows  the  shape  of  the  applied  force 
for  each  impulse,  (c)  shows  that  the  two  total  impulses  are  in  fact  different  by 
about  a  factor  of  two. 
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Figure  12  The  TRIM  process  can  be  used  to  resolve  discrete  forces  applied  to  the  NIBS,  (a) 
shows  the  deflection  trace  corresponding  to  (b)  the  time-varying  force  traces  for 
2  discrete  impulses  of  IVOmN  for  1  sec  each,  (c)  shows  the  resulting  total 
impulse  as  a  function  of  time. 

Figure  13  Plots  of  (a)  the  applied  and  derived  forces  as  a  function  of  time  for  an  arbitrary 
impulse,  and  (b)  the  corresponding  applied  and  derived  total  impulses. 

Figure  14  Comparison  of  the  analytical  model  with  the  experimental  measurements  for 
multiple  pulses  of  varying  pulsewidth.  The  plots  show  the  sequence  of  steps  in 
the  TRIM  process,  (a)  The  balance’s  deflection,  (b)  velocity,  and  (c)  acceleration, 
(d)  The  time  varying  force,  and  (e)  the  total  impulse  as  a  function  of  time. 
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Table  1 


Parameter 

Value 

Error  (+/-) 

Moment  of  Inertia,  I 

0.03184  kgm" 

0.00010  kg  m^ 

Damping  Coeffieient,  C 

0.00023  kg  m"/s 

0.00002  kg  m^/s 

Spring  Constant,  K 

0.20034  N/m 

0.00150  N/m 

Period,  T 

2.49  s 

0.01  s 

Angular  Frequeney,  ODeff 

2.523  s'" 

0.010  s'" 

Damping  Ratio,  5 

0.0090 

0.0009 

Moment  Arm,  r 

0.4635  m 

0.0005  m 
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Figure  1 
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Figure  2 
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